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We develop a quantum field theoretical framework to analytically study the three-body con- 
strained Bose-Hubbard model beyond mean field and non-interacting spin wave approximations. It 
is based on an exact mapping of the constrained model to a theory with two coupled bosonic degrees 
of freedom with polynomial interactions, which have a natural interpretation as single particles and 
two-particle states. The procedure can be seen as a proper quantization of the Gutzwiller mean 
field theory. The theory is conveniently evaluated in the framework of the quantum effective ac- 
tion, for which the usual symmetry principles are now supplemented with a "constraint principle" 
operative on short distances. We test the theory via investigation of scattering properties of few 
particles in the limit of vanishing density, and we address the complementary problem in the limit 
of maximum filling, where the low lying excitations are holes and di-holes on top of the constraint 
induced insulator. This is the first of a sequence of two papers. The application of the formalism to 
the many-body problem, which can be realized with atoms in optical lattices with strong three-body 
loss, is performed in a related work |14J. 
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I. INTRODUCTION 

Lattice theories with constrained bosons have proven 
to be a powerful description of various spin models and 
strongly correlated systems in condensed matter physics 
[T]. On the other hand, such theories with constrained 
lattice bosons have recently arisen naturally in effective 
models for experiments with cold atoms in optical lat- 
tices. In the presence of large two-body and three-body 
loss processes, bosons in an optical lattice are described 
on short timescales by a model with two-body and three- 
body constraints respectively \%\ i3j . The behavior of the 
Bose gas is changed drastically; for example, in the case 
of the three-body constraint, the creation of an attrac- 
tive Bose gas with atomic (ASF) and dimer superfluid 
(DSF) phases is possible While the possibility of such 
an ASF-DSF transition has been predicted earlier in the 
context of continuum attractive Bose gases near Feshbach 
resonances [U E] , the constrained lattice system offers an 
intrinsic stabilization mechanism to observe such a phe- 
nomenology in experiments. This serves as one motiva- 
tion to study such models theoretically in more detail, in 
particular exploring the consequences of the presence of 
the constraint. We also note that by the same dissipative 
blockade mechanism, constrained models with fermions 
may be created [SI [7] . 

Here our goal is to describe the physics of a constrained 
boson system beyond a mean-field plus spin wave ap- 
proach (see e.g. [T]). In order to do this, we find an ex- 
act mapping of the original constrained bosonic Hubbard 
model to a theory of two coupled unconstrained bosonic 
degrees of freedom which interact polynomially. The re- 
sulting theory is conveniently analyzed in the framework 
of the quantum effective action, which makes it possible 
to study both thermodynamical and dynamical proper- 



ties of the system via various many-body techniques. As 
a consequence, we can demonstrate several remarkable 
features of the three-body constrained attractive Bose 
lattice gas, which are uniquely tied to this constraint and 
not treated properly within a simple mean field plus spin 
wave approach. In particular, we show the emergence of 
an Ising quantum critical point on the phase transition 
line between atomic to dimer superfluid phases, which 
generically is preempted by the Coleman- Weinberg mech- 
anism j8j where quantum fluctuations drive the phase 
transition first order, rendering the correlation length fi- 
nite [U El IMTT) . We also show the presence of a bicritical 
point [12] in the strongly correlated regime at unit filling 
of bosons, which is characterized by energetically degen- 
erate orders. In our case this corresponds to the coex- 
istence of superfluidity and a charge density wave. Fur- 
thermore, quantitative effects of quantum fluctuations on 
the position of the phase boundary can be investigated 
systematically. 

The formalism we develop here has much broader ap- 
plications than the bosonic lattice gas with a three-body 
constraint. In particular, it could be used to treat sys- 
tems with effective constraints arising from large interac- 
tion parameters, and it is also applicable to constrained 
fermionic models, which could arise due to strong loss 
features in 3-componcnt fermion gases [7j [13] . As a re- 
sult, we focus in this paper on presenting the quantum 
field theoretical construction in detail, and give bench- 
mark calculations for this method. Its application to the 
many-body attractive lattice Bose gas with three-body 
constraint is left to [14., where we discuss in detail the 
results outlined above. The results of the present and the 
related paper |14j are summarized in Il5j. where we also 
indicate how to probe our findings in experiments with 
ultracold atoms. 



2 



This paper is organized as follows. We begin by re- 
viewing the microscopic derivation of the three-body con- 
strained model via a dissipative blockade mechanism in 
Sec. [n] In Sec. |III| we give an overview of our formalism 
and a comparison to existing methods for treating con- 
strained models. Sec. IIVI contains the central result of 
this paper, the mapping of the constrained model to a 
coupled boson theory already anticipated above. In Sec. 
[V] we apply the formalism to the two limits where true 
many-body effects are absent, n = and maximum filling 
n = 2. There our boson model reduces to Feshbach-type 
models, which we analyze in terms of Dyson-Schwinger 
equations [16j . We perform the above mentioned bench- 
mark calculations at n = 0. Furthermore, we investi- 
gate hole/di-hole scattering and bound state formation 
at n = 2, and present a fourth order perturbative cal- 
culation for the dimer-dimer interaction. This will be 
needed as an input for the many-body theory in |14j . 
Our conclusions are drawn in Sec. IVII 



II. DERIVATION OF THE CONSTRAINED 
MICROSCOPIC MODEL 

In this section, we review how the three-body hardcore 
constraint for the Bose-Hubbard model emerges from 
strong three-body loss [3] . 

We consider bosons in the lowest Bloch band of an 
optical lattice, which are described by the Bose-Hubbard 
Hamiltonian {h= 1), 

where ai(aj) is the bosonic annihilation (creation) opera- 
tor at site i. J is the hopping rate, n the chemical poten- 
tial, and U the onsite interaction. The convention 
first sums over all sites i, and then over the neighbour- 
hood of each i spanned by sites j, J2{i,j) = J2iJ2{j\i)- 
This model is valid in the limit where J, Un <C uj, where 
w is the separation between Bloch bands and n is the 
mean density. 

Three-body loss in this system is due to inelastic colli- 
sions of three atoms, two of which form a deeply bound 
molecule. Together with the third atom, molecule for- 
mation is compatible with energy and momentum con- 
servation, unlike the case of two-particle collisions. Since 
the binding energy of the molecule typically strongly ex- 
ceeds optical lattice depths, the resulting kinetic energy 
of the products couples them to the continuum of un- 
bound states, thus leading to their escape from the lat- 
tice. This picture allows us to write down a zero tem- 
perature master equation in the Markov approximation, 
which in the simplest approximation neglects loss arising 
from particles on neighbouring lattice sites. This master 



equation is given by 

P = -i {HoffP - pHl^) + ^ E 2«'P«r' (2) 

i 

TJ TT \ " t 3 3 

i 

where the decay rate 73 can be roughly estimated from 
the experimentally measured continuum loss rate via the 
usual Wannier construction [3]. We have absorbed the 
Schrodinger type terms of the dissipative evolution into 
an effective Hamiltonian with imaginary, and thus decay, 
term. The remaining recyling term couples sectors in 
the density matrix with particle number TV, TV — 3, ... and 
ensures a norm conserving time evolution. 

We are now interested in the limit of strong loss, 
73 3> J, t/, which suggests a perturbative expansion in 
1/73. The scale 73 only couples states with three and 
more particles per site. We therefore define a projector 
P onto the subspace with at most 2 atoms per site, and 
its complement as Q = 1 — P. To second order in pertur- 
bation theory we then find the projected Hamiltonian 

i?P,eff « PHbkP + PHBKQiQHBKQy^QHBHP 

= PHBHP-JJ2P4c^P, (3) 

i 

and the corresponding Master Equation is 
pp = -i {p[p,cspp - ppHlcff) + r E Pc^pp4 ""P, (4) 

i 

with Pp ~ PpP. The effective decay rate and jump 
operators are given by 

j2 

F « 12 — , c^a^E 0^/72. (5) 

The respective terms in -ffp.cff have simple interpreta- 
tions: The leading term describes the coherent dynamics 
of lattice bosons, but with the constraint of not popu- 
lating a single site with more than two particles. This 
projected part H = PHbhP can thus be written as a 
3-body constrained Bose-Hubbard Hamiltonian 

H = -JY^ - pYI ""J"' + 51 "J 

(ij) i i 

«r = 0. (6) 

The leading correction is imaginary and describes parti- 
cle number loss. The decay rate F I/73 is, however, 
very small in the considered limit. Consequently, over 
timescales r = 1/F, one realizes indeed the physics of the 
Bose-Hubbard model with 3-body hardcore constraint. 
In 0], we have shown that e.g. in atomic Cesium systems 
close to the zero crossing of the scattering length, the loss 
rate 73 is the dominant energy scale. There, we have also 
analyzed the many-body dissipative dynamics of Eq. ([2]) 
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in the regime described by Eq. ^ with exact DMRG 
methods in one spatial dimension, including the specifi- 
cation of a scheme with which the ground state of PHP 
can be reached. This motivates a more detailed analyti- 
cal investigation of the zero temperature phase diagram 
of Eq. As mentioned in the introduction, such a sce- 
nario with dominant three-body loss also arises naturally 
in three-component fermion systems close to a Feshbach 
resonances |13j . where it arises due to the proximity to 
rapidly decaying Efimov state. 



III. OVERVIEW 

In this section, we provide an overview and discussion 
of our construction and compare it with existing theoret- 
ical approaches. 

The starting point is a truncation of the onsite bosonic 
Hilbert space to three states corresponding to zero, sin- 
gle and double occupancy. Following Altman and Auer- 
bach [24j, we introduce three operators creating these 
states out of an auxiliary vacuum state. The operators 
are not independent but obey a holonomic constraint, 
such that one of the degrees of freedom can be elim- 
inated. We propose a resolution of this constraint in 
a fashion that (i) produces polynomial interactions be- 
tween the two remaining operators and (ii) allows us to 
interpret them as standard bosonic degrees of freedom. 
The Hamiltonian written in terms of these bosonic oper- 
ators is an involution on the physical Hilbert space (as- 
sociated to onsite occupation < 2), i.e. all matrix el- 
ements coupling into the unphysical sector with onsite 
occupation > 2 vanish. Consequently the partition sum 
for this Hamiltonian decomposes into a physical and an 
unphysical part. In order to distinguish the respective 
contributions, we choose to calculate the quantum effec- 
tive action from a functional integral representation of 
the partition sum. This is the most general polynomial 
in the operators of the theory which is compatible with 
the symmetries of the microscopic Hamiltonian, respec- 
tively the microscopic action. We argue that the general 
restrictions imposed by symmetry are supplemented by a 
further new "constraint principle" , the requirement that 
the contributions to this polynomial be compatible with 
the constraints imposed by the microscopic Hamiltonian 
- the evaluation can proceed as in a standard polyno- 
mial boson theory. This opens up the powerful toolbox 
of modern quantum field theoretical methods for calcu- 
lations in onsite constrained models, and the effective 
action provides a unified framework to treat the strong 
correlations at short distances characteristic to lattice su- 
perfluids as well as the deep infrared sector of the theory, 
where questions on the nature of the phase transition can 
be addressed. Similar to symmetries, the restrictions im- 
posed by the constraint on the full theory leverage over 
from the microscopic theory. However, unlike symme- 
tries, the relevance of the constraint depends on scale, 
being restrictive on short distances, while on long dis- 



tances power counting arguments lead to an effectively 
unconstrained though interacting spin wave theory. 

Our approach may be conceived as a proper quanti- 
zation of the Gutzwiller mean field theory: The con- 
struction implies that the zero order contribution to the 
thermodynamic effective potential is the Gutzwiller mean 
field energy. The quadratic fluctuations reproduce spin 
wave theory. However, due to the exact nature of the 
mapping we can also assess the effects of interactions sys- 
tematically. At the same time, this shows that while the 
mapping is exact, it builds on the knowledge of qualita- 
tive properties of the ground state of the system, such as 
the symmetry breaking pattern. In other words, given 
this pattern, we can systematically construct the excita- 
tions and their interactions on top of it. The bias intro- 
duced by the assumptions on the ground state is actually 
shared with all analytical approaches to many-body sys- 
tems involving symmetry breaking. 

One established way to deal with constrained bosonic 
theories with a finite number of onsite states would be 
to map it to the corresponding spin model (spin 1 in 
our case), and subsequently analyze this model. In prac- 
tice, however, the analytical evaluation of the spin model 
beyond a free spin wave level is technically challenging. 
Furthermore, the spin language is not particularly advan- 
tageous and physically intuitive for models which origi- 
nate from a constrained boson theory, cf. Sect |IVD| We 
see a central advantage of our formalism vs. such a spin 
model approach in the straightforward interpretation of 
the two bosonic fields in physical terms, which are close 
to the original degrees of freedom. For example, at low 
density n > and close to maximal filling n < 2 they 
may be interpreted as atoms and bound states of these, 
dimers, resp. holes and di-holes, which are indeed ex- 
pected to form the dominant low energy excitations in 
the respective situations. At intermediate densities hy- 
bridization takes place. The microscopic Hamiltonian 
reduces to Feshbach-type models, with the characteristic 
splitting term of diniers/di- holes into atoms. The role of 
the detuning is now played by the interaction term U . We 
may view this result as a built-in Hubbard-Stratonovich 
transformation on the level of the Hamiltonian, which 
importantly respects the constraint. At this point, we 
note an important difference of our approach to the one 
by Dyson and Maleev [T7] , which map a spin model to a 
model with a single bosonic degree of freedom. The gen- 
eralization of a Hubbard-Stratonovich transformation to 
such a model would be problematic. 

Other related approaches to constrained models have 
been put forward in the context of slave boson theories for 
the description of strongly correlated fermions [HI [12] , or 
in the Schwinger-boson or Holstein-Primakoff approach 
to spin models . While these formulations are in princi- 
ple exact, the practical implementation of the constraint 
is done on an approximate level only |19j , or in a fashion 
that makes it difficult to assess coupling constants of an 
effective low energy theory or thermodynamic quantities 
[18) . Our approach is tailored to make these accessible. 
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Furthermore, we stress that our exact implementation of 
the constraint differs from the conventional approximate 
treatment via expansion of the square root as done e.g. 
for the case of large spin S in the Schwinger-boson or 
Holstein-Primakoff approaches [Ij. We provide concrete 
benchmarks for our procedure through the analysis of the 
vacuum problem of a few scattering particles. While our 
approach yields the correct non-perturbative Schrodinger 
equation for two-particle scattering and the right coeffi- 
cient for an induced nearest-neighbour dimer-dimer inter- 
action, the square root expansion would produce wrong 
prefactors. 

Clearly, going beyond the mean field approach in 
constrained bosonic models is nowadays possible mak- 
ing use of powerful numerical techniques such as Quan- 
tum Monte Carlo [^DH2^ or variational simulations [53]. 
However, often a more analytical understanding of the 
physics leading to a certain quantitative effect is desir- 
able. Furthermore, as we have mentioned above, our for- 
malism can be extended to fermion systems, which is not 
straightforward with numerically exact techniques. 



IV. MAPPING THE CONSTRAINED MODEL 
TO AN INTERACTING BOSON THEORY 

In this section we derive the mapping of the bosonic 
model with three-body hardcore constraint Ec^. (|6| to 
an unconstrained model with two bosonic degrees of free- 
dom, coupled via polynomial interactions. This construc- 
tion depends on the filling of the lattice. We concentrate 
on the "vacuum limit" of zero density and temperature 
first, where one deals with a few particles and many-body 
effects, such as spontaneous symmetry breaking, are ab- 
sent. The generalization to arbitrary filling < n < 2 is 
performed in Sec. |IVC[ 



A. Reduction to three on-site states 

We start from the 3-body constrained Bose-Hubbard 
Hamiltonian Eq. The constraint makes it possible 

to restrict the Hilbert space on each site to three states 
with occupation 0,1,2. We introduce three operators that 
create "particles" in these states. Such a slave-boson type 
procedure has been proposed previously by Altman et al. 
[24] in the context of the repulsive Bose-Hubbard model, 
where it constitutes an approximation. Here this step is 
exact: 
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I") = ^a.Jvac) 



1 



a] "I vac) 



0,1,2. (7) 



The operators ta^i are so far only defined by their action 
of the "vacuum" state |vac) and we discuss their commu- 
ation relations below. They obey a holonomic constraint. 
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a=0 



(8) 



In the space spanned by |0), |1), |2) we may express the 
creation operator as 



(9) 



Hence we can express the original Hamiltonian equiv- 
alently in terms of the new operators, defining non- 
hermitian kinetic operators 



(21) _ 



+ 1 / 



(10) 



and hermitian potential energy operators 

"l.i ~ "2,1 ~ ^2.1*2,4, (11) 



we obtain 



(12) 



+^/2(xfl)iff°)t+^(lo)^(2l)t)]^ 
-ffpot = -fJ-^ni + U^n2,i, fii ^ fii^^ + 2n2,i 



We observe that the local terms become simple 
(quadratic) in this representation, while the nonlocal 
terms are fourth order in the operators, giving rise to 
"kinematic" interactions. The onsite interaction part 
therefore can be treated exactly, while the complexity of 
the problem is now encoded in the hopping term. This is 
reminiscent of the conventional Gutzwiller (mean field) 
approach to the Bose-Hubbard model. Indeed, imple- 
menting a Gross-Pitaevski-type mean field theory for the 
above Hamiltonian by formally replacing the operators 
by complex valued amplitudes ta,i — > /a,i, the above 
Hamiltonian operator reduces to the Gutzwiller energy 
expression i?GW = {'^\H\'ip) for the case that the wave 
function — is truncated to the three low- 

est Fock states on each site, — X]q=o 
this case, the holonomic constraint reduces to the nor- 
malization condition for the wave function, i{tp\^)i — 
HaPaJcS = IVi. 

It may be tempting to try to develop a many-body 
theory including the description of spontaneous symme- 
try breaking at finite density for the above Hamiltonian 
directly in terms of the ta operators in the sense of a 
Bogoliubov-type theory on top of the Gross-Pitaevski 
mean field, via a replacement of the type t^ i = fa,i+Sta^i 
familiar from low density continuum theories. Such a 
procedure, however, leads to severe consistency problems 
when encompassing the full range of densities allowed by 
the three-body hardcore constraint, < n < 2. There- 
fore, we first focus on the "vacuum limit" n = where 
no spontaneous symmetry breaking is present, and where 
the Hamiltonian ( [l2| describes the physics of a few scat- 
tering particles. The generalization to arbitrary density 



is performed in Sec. IV C 
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B. Implementation of the Holonomic Constraint: 
Interacting Boson Theory 

The holonomic constraint is now used to eliminate one 
of the operators. In the limit n = 0, all the amplitude 
resides in the zero-fold occupied state, and the mean field 
vacuum is described by = Hi il'^^'^)- thus ehm- 
inate the operators io,i. The remaining two operators 
describe excitations on top of this mean field vacuum, 
and will have a natural interpretation in terms of atoms 
and dimers. At this point, it is transparent that our 
construction builds on the proper choice of the qualita- 
tive features of the physical vacuum. This prerequisites 
a certain understanding of the physics, and introduces 
a bias in our construction. The subsequent construc- 
tion is however exact, and can be used to quantitatively 
calculate properties of the system. The situation is actu- 
ally similar to the treatment of interacting Bose gases in 
the continuum, where a condensation of the particles in 
the zero mode is assumed to expand in the fluctuations 
around this mean field. 

The holonomic constraint can, in principle, be imple- 
mented by the replacements 



4,.e''^Vl - "1,^ - ^2.^, (13) 



The holonomic constraint is local and only restricts 
the amplitude of fg.ij therefore we introduce a phase- 
amplitude representation |<o,i| exp 1^?^. We first discuss 
the role of this phase and then turn to the more interest- 
ing question of the amplitude. 

Inserting the replacement into Eq. (12), we observe 
that a local redefinition of ti and t2 as 



t^ A ^ t-i i — t 



(14) 



for all i completely removes the phase from the Hamilto- 
nian. Therefore, we can work in this rotated frame from 
the outset, and simply consider the two complex valued 
operators ?i;2,i. We will drop the tilde in the following. 

Now we study the amplitude. Obviously the square 
roots are impracticable for any field theory calculation 
where one has to work with polynomials in the field op- 
erators. However, on our subspace, the matrix elements 
of both (1 — nij — n2j)^^^ and (1 — nij — n2,j) are the 
same: either 1 or 0. Consequently, on the subspace we 
may replace 



X,, 



1 - ni- 



- ^^2,^• 



Xiti . 



(15) 



Note that the second expression could also be replaced 
by (1 — ni i — n2,i)^^^ii,i — >■ ti i, since the action of this 
operator is nonzero only if originally there is an atom on 
the i-th site, but when this atom is annihilated by ti^i 
with an empty state left, the action of the square root is 
simply unity. However, for hermeticity issues we prefer 



Formally, we may justify the replacement of the square 
root by a polynomial by the following formula: Consider 
a linear operator X with the property X"^ = X. Then for 
a function / of this operator one has, using the Taylor 
representation. 



fix) 



= f{0){l-X)+X 



= f{0){l-X) + .f{l)X. 



(16) 



In our case, we have X — 1 — = ni^i+fi2,i for all i, and 
f{X) = \/l ~ X. Indeed, X'^ = X. Seen as a function of 
X,^l-X,we have /(X,) ^ VX, ^ /(0)X, + /(1)(1 - 
Xi). Since also Xf — Xi the latter result would have 
been obtained from the Taylor representation of f{Xi) 
around 0. The auxiliary operator X is introduced to 
circumvent an expansion of the square root around 0, 
but leads to the same result. No approximation has been 
used here. 

Having implemented the constraint, we are now going 
to show that the remaining operators ti , t2 can be treated 
as standard bosonic operators. Consequently the Hamil- 
tonian (or the corresponding action) with the above re- 
placements will lend itself for a treatment with well es- 
tablished field theoretic methods. To show that we may 
interpret ti, t2 as bosonic operators, we assume a bosonic 
Hilbert space for atoms ti and dimers t2 at each site i, 
which reads T-Li — {\ni) Ito^)}, Ui^rui = 0, 1, .... The com- 
plete Hilbert space is "H = Hi^i- We divide the onsite 
Hilbert spaces into a physical subspace Vi and an unphys- 
ical one Ui, Hi — Vi © Ui ; the subspaces are orthogonal 
by construction. The physical subspace is spanned by 
the combinations 



7'z = {|0,)|0,),|l,)|0,),|0,)|l,)}. 



(17) 



to work with the variant in Eq. ( 15 1 



In our construction of introducing atom and dimer oper- 
ators, the state with two atoms on one site is represented 
as one dimer |Oi)|li); |2i)|0i) instead is already part of 
U,. 

A first observation is that standard bosonic operators 
have the same action on the physical subspace as the orig- 
inal operators ti, t2 defined via Eq. (0, since the bosonic 
•y/n-enhancement {W\n) = y/n + l\n -I- 1) for bosonic op- 
erators h^) is either or 1 on the physical subspace. The 
assumption of ti , t2 being bosons is thus consistent on 
the physical subspace. 

with the 
le goal to 



Next we consider the Hamiltonian Eq. ( 12 
constraint implemented via Eq. (151, with t 



show that the time evolution generated by this Hamilto- 
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nian does not couple the two subspaces. It reads 



n2.i~f^ni,i, (18) 



In the kinetic term, the first expression describes the con- 
ditional hopping of single atoms, the second represents 
the exchange of a dimer and an atom on neighbouring 
sites, and the last one describes the conditional bilocal 
splitting and recombination of a dimer into atoms. It 
can be easily shown that H maps physical on physical 
states, and unphysical on unphysical ones, while there 
are no transitions between the subspaces generated by 
(18). For that purpose it is sufficient to check that 



Hi^.j\ni)\mi)\nj)\m.j 



(19) 



is in ViVj if the initial state is in ViVj (a simple 9 dimen- 
sional space). Since this means that all matrix elements 
{u\Hij\p) = 0, this is sufficient to conclude that starting 
in U the mapping will be into U, since for the hermitian 



H one has (p|i?r 



\Hij\p)* ~ 0. In consequence 



H can be written in the form 



H = Hp®li 



(20) 



In other words, "P = Hi ^ — Yli^i invariants 
under application of H, and thus also repeated applica- 
tion does not lead out of the subspaces. Therefore, we 
also have 



expi-^H) = exp(-/?i7p) (g) ![/ -f Ip exp{~(3Hi 
The partition sum is the given by 



(21) 



Z = Trexp(-/3i7) 



(22) 

\P)' 





exp{-l3Hi, 

= E(p| eM'PHp)\p) + eM-PHu)\u). 
{p} {«} 

Thus we get contributions from both the physical and 
the unphysical part of the Hilbert space. However, the 
key point is that they do not mix; thus the answers found 
for the physical part will be correct, and we only need to 
find the criterion to discriminate the physical from the 
unphysical part of the partition sum. 

This issue is addressed in the last step of the construc- 
tion. Indeed, such a setting is provided by using the 
effective action to encode the physical information of the 
theory, as we will now outline. First we represent the par- 
tition function as a Euclidean functional integral, which 



is straightforward as we simply have to quantize a theory 
with two coupled bosonic degrees of freedom [25] : 



Z = J VtiVt2CXp-S[ti,t2], 



(23) 



S = 



dT 



drti 



^2 i^Tt2A 



H[tiM 



where H is the Hamiltonian above, however to be inter- 
preted in the Heisenberg picture with (imaginary) time 
dependent fields, and the fields are now classical fluctu- 
ating variables. S is the classical Euclidean action. In 
the next step we introduce a source term in the partition 
function. 



Z[ji,j2] = j Vt{Dt2exp-S[t^M 



(24) 



/ d'T^{il,itl,i + il,it2,i + C.C), 

i 

Z[ji=j2=0]. 



32 

The source terms introduce linear terms in ti^2, which 
mix the physical and the unphysical sectors. Since ji_2 
is only used in a pivotal sense to generate the correla- 
tion functions upon functional differentiation, and set to 
zero at the end of the calculation, this does not pose any 
conceptual problems. The situation is analogous to the 
effect of the source term on the symmetries of the theory, 
which are broken explicitly for nonzero sources. 

The effective action is defined as the Legendre trans- 
form of the free energy W[j] = logZ[j] (we introduce the 

shorthands x = (ti, 4, t2, 4)> i = {h-, i\, 32, il) [26,: 



■T 

] X, X 



6j 



(25) 



where x — ix) is the field expectation value or the "clas- 
sical" field. By the Legendre transform, the active vari- 
able is changed from j to x- The effective action has the 
following representation in terms of a functional integral. 



exp-r[x] = J VSxcxp-S[x + Sx]+ J 



j 5x, 



J = 



Sx 



(26) 



where Sx = X ^ X- The last identity is the full quan- 
tum equation of motion, and the equilibrium situation we 
are interested in is specified hy j = where no mixing 
between the physical and the unphysical sector occurs. 
When fluctuations are unimportant, the integration over 
the Sx can be dropped and the above equation reduces 
to r[x] = S[x], i.e. the quantum effective action reduces 
to the classical one. 

The effective action expresses the theory in terms of the 
fields X- The vertex expansion generates the one-particle 
irreducible (IPI) correlation functions, 

r[x]==E7i / ^il....,^li^^^■■■'^l)x^Axl) ■ ... ■ Xh{xi)- 



(27) 
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Usually, the coupling coefficients of the expansion are 
only restricted by the symmetries of the theory - the ef- 
fective action is the most general polynomial in the fields 
X which is compatible with the latter. Thus, formulating 
the theory in terms of physical objects - the fields x ~ 
offers the advantage of directly leveraging the power of 
symmetry considerations from the microscopic (or clas- 
sical) to the full quantum level. 

In complete analogy, we can make use of the restric- 
tions present in the microscopic Hamiltonian when com- 
puting the quantum effective action. Since, as we have 
shown above, no couplings mapping from U ^ V are gen- 
erated, we may write down the most general form for the 
effective action for the physical sector of the theory by 
directly excluding couplings which would violate this con- 
straint. In practice, this concerns processes which change 
the on-site occupation number. For example, a process 
which involves creation of a dimer on site i must be ac- 
companied by an appropriate constraint that the site be 
empty prior to the process. Thus, the operator t\ must 
always appear in the combination t^Xi. Furthermore re- 
quiring hermeticity of the terms appearing in the effective 
action, we conclude that the effective nearest-neighbour 
dimer hopping term is of the form J^st^ ^XiXjt2j- In 
the practical calculation, we may restrict ourselves to 
the computation of the coefficient which is simplest to 
extract - obviously the quadratic one. 

In sum, we have obtained the following result: the 
usual symmetry constraints on the quantum effective 
action are now supplemented by a further fundamen- 
tal principle, namely the restrictions present in the mi- 
croscopic theory which originate from the hardcore con- 
straint. 



C. Arbitrary density 

The construction presented in the last section focused 
on the zero density limit, describing the scattering of few 
particles in the absence of many-body effects. At a finite 
density, the low temperature physics of bosons is char- 
acterized by the spontaneous breaking of global phase 
rotation symmetry U{1). The ground state exhibits a 
condensate mean field, which has to be incorporated in 
the theoretical description of the system. One customary 
approach is to quantize a theory with degrees of freedom 
h = s + 5b via the path integral, where 5b is the fluctu- 
ation around the classical field s. However, in our case 
this procedure does not work, since possible values of the 
mean field lie on the compact interval [0, 2], and due to 
the nonlinear nature of the constraint. To this end we 
follow Huber et al. '27j , who implement the procedure on 
the mean field plus spin wave level. We will see that our 
treatment of the constraint can be applied also in this 
case, such that we arrive at an exact formulation of the 
problem in the presence of spontaneous symmetry break- 
ing. The procedure consists in first introducing the mean 
field via a unitary rotation in the space of operators (Eq. 




FIG. 1: Rotation to a new ground state: The vacuum 
states n = 0(2) are described by the red mean field vec- 
tor pointing in positive (negative) z direction. The fluctu- 
ations (black arrows) form a coordinate system in which the 
direction coUinear to the mean field vector is eliminated via 
the implementation of the constraint. All mean field vec- 
tors not in the z direction describe a homogeneous superfiuid 
ground state. A situation analogous to the one for the vacua 
is achieved via an appropriate rotation of the coordinates for 
the fluctuations. 



( 29 1 below) , and then quantizing this theory of operators 
which are free of expectation values via the path integral. 
In this way, we can obtain a picture which is fully con- 
sistent with general features of many-body theories with 
[/(I) symmetry, in particular, we can derive Goldstone's 
theorem within our framework |14) . 

Here we concentrate on homogeneous ground states, 
noting that the implementation of spatially dependent 
order parameters - such as a charge density wave - is 
straightforward. One then simply has to use rotation 
matrices which vary from site to site. Such a situation 
will be encountered in [H]. 

Our treatment of the vacuum problem at n = started 
from the idea that in this case, all the amplitude resides in 
the zero-fold occupied state, and that fiuctuations around 
this state have to be considered. The excitations on top of 
this "mean field" vacuum, defined as \Vt) = Ili^ il'^^'^)' 
then turned out to be single atoms and dimers, respec- 
tively, as expected intuitively, and we have formulated 
the corresponding quantum field theory to describe their 
scattering properties. 

In the many-body problem, we proceed in complete 
analogy by first introducing a mean field vacuum. A 
general homogeneous mean field vacuum may be written 
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as 



= n(E^"^'^p(*"'^)i")o (28) 

i Q 

i a i 

The introduction of 6j ^ as the new vacuum creation op- 
erator implies the need for a redefinition of the remain- 
ing two degrees of freedom. Such a transformation is 
performed via a two-parameter unitary rotation, whose 
rotation angles are chosen such that the new operators 
fluctuate around the new vacuum state and do not fea- 
ture expectation values (cf. Fig. [T]), 

= (^fl^x)a/3^L (29) 

with the explicit form of the rotation matrices 

cos6'/2 sine'/2e2'"^ 
Rb = \ 1 I , (30) 

-sin6'/2e-2i<^' o cos 61/2 

1 

i?x = I cosx/2 -sinx/2e''^ 
sinx/2e^''^ cosx/2 

A finite 0{x) corresponds to a finite amplitude in |2)(|1)). 
The precise relation is 



ro ~ cos 0/2, ri = sin 0/2 sin x/2, r2 = sin 0/2 cosx/2. 



(31) 



The strategy is to first rotate to the new mean field state 
by inverting the unitary matrix in Eq. ( 29 ) and sub- 
sequently implement the constraint. Analogous to the 
procedure in the physical vacuum, we may now eliminate 
the operator b^ which is chosen to include the expectation 
value. We note that the local rotating frame transforma- 
tion ( 14 1 can be applied also here, showing that the phase 
of 5o is irrelevant for the Hamiltonian. Consequently we 
can implement the constraint by the formal replacement 



blfi^^.^X,. (32) 



The second expression is simply a rearrangement of the 
holonomic constraint. The resulting bosonic Hamilto- 
nian, which is then quantized by means of a functional 
integral, is rather complex, and its explicit form and anal- 
ysis are discussed in |14j . However, it exhibits a simple 
structure. 



i?GW is the Gutzwiller mean field energy and -ffsw de- 
scribes the quadratic spin wave theory ^ . The corrections 
to the mean field phase diagram, as well as nontrivial ef- 
fects in the deep infrared physics which we analyze in 
[T3] are not captured at this quadratic level. They are 
all encoded in the interaction part i?int- 

At this point, let us compare our findings to the work 
of Huber et al. [27^ • We have verified explicitly that the 
quadratic part of the Hamiltonian coincides with the spin 
wave or Bogoliubov theory obtained in that work, though 
the authors use a different prescription for the resolution 
of the constraint for a single 6o operator. The reason is 
that whenever the operator 5o appears in multiplicative 
combination with another operator 6i , 62 , the replace- 
ment prescribed by the first expression of Eq. ( 32 1 gives 



rise to at least cubic terms neglected in the spin wave ap- 
proximation. In contrast, when the combination in the 



GW 



sw 



(33) 



second expression of Eq. ( 32 1 appears, one simply has to 
rearrange the constraint and no difference between the 
two approaches appears. 

The many-body problem is completely specified only 
upon indicating conditions which determine the two ro- 
tation angles 0, x- In [Mj we show that they are fixed by 
gap equations which emerge as a consequence of Gold- 
stone's theorem. Furthermore, the chemical potential 
is fixed via the equation of state when calculations at 
fixed filling are intended. These conditions are exact but 
implicit, and can be resolved approximately only. The 
simplest approximation reproduces the classical or mean 
field solution as anticipated above, but it is possible to 
go beyond this simple scheme with our setting |14| . 



D. Relation to a Spin-1 Model 

From the availability of only three onsite states, it is 
clear that there exists a mapping of the 3-body con- 
strained Bose- Hubbard model to a spin-1 model. Nev- 
ertheless, it seems advantageous to us to work with the 
above construction of mapping to a coupled boson theory. 
The reason is twofold. 

First, the practical analytical analysis of spin models 
beyond the standard mean field plus non-interacting spin 
wave approximation is technically very hard. Indeed, 
usually one resorts to introducing an artificial smallness 
parameter (inverse number of field components 1/N , in- 
verse total spin 1/S), such that a Gaussian theory be- 
comes exact in the limit N,S — >■ 00 , and organizes an 
expansion scheme around this noninteracting fixed point 
(1/A^, 1/5' expansions). This is then followed by a con- 
tinuation of the results to the physical system of inter- 
est, where N and S are typically small. However, in 



The phase relation between the amplitudes fa = Vae^^" ,0a = ^ A linear contribution, as naively expected in the expansion about 
aij> is a consequence of spontaneous symmetry breaking in Fock the condensate, does not occur due to Goldstone's theorem, see 

space. is the spontaneously chosen phase of the condensate. Ref. I14| . 
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our model 5 = 1, and most of the effects which we dis- 
cuss here and in |14j - ranging from the non-perturbative 
formation of the dimer bound state in vacuum over cor- 
rections to the phase boundary to the true nature of the 
phase transition - are not accessible to leading order in 
the abovementioned schemes, since they are all rooted in 
the intrinsic non-linearities of the theory. 

Second, typically the direct mapping of a bosonic the- 
ory with hardcore constraint yields a rather complicated 
effective microscopic spin Hamiltonian, which further 
complicates an analysis in terms of spin degrees of free- 
dom. For example, in our case the corresponding spin 
model would feature cubic and quartic bilocal spin in- 
teractions ~ sfsfsj + h.c.,^ sfsfsJSj + h.c. with in- 
teraction constants of the same order as the quadratic 
terms [27] . These terms break the rotation symmetries 
typically present in generic Heisenberg models for mag- 
nets. Physically, the appearance of such terms has to 
be expected, since none of these rotation symmetries are 
present in the hardcore boson model. Clearly, mapping 
the constrained model to a theory of unconstrained cou- 
pled boson degrees of freedom, which find a natural in- 
terpretation in terms of single particle and bound state 
degrees of freedom, is closer to the physics of the boson 
theory with hardcore constraint. 

Reversely however, we emphasize that the class of spin 
models which can be readily encompassed within our for- 
malism is made up of Heisenberg XX models in external 
fields, with possible extensions to XXZ models with small 
anisotropy. The power of conventional field theory tech- 
niques can thus be applied to such models. 

V. FLUCTUATIONS IN THE VACUUM 
PROBLEM 



play a role. However, a mathematically wrong or only ap- 
proximate implementation of the constraint will produce 
a wrong scattering equation as we see explicitly. Second, 
we use the language of Feynman diagrams to explicitly 
calculate the two-dimer interaction strength up to fourth 
order in the perturbative regime J/\U\ ^ 1. This pro- 
vides another benchmark for our formalism. Third, we 
consider hole scattering in the limit n = 2 and discuss 
the formation of a di-hole bound state, which exhibits 
properties different from the dimer bound state at n = 0. 

Below we introduce the formalism used to do concrete 
calculations, the Dyson-Schwinger equations, demon- 
strating how to use the powerful methods of quantum 
field theory in our problem. These equations provide 
an exact hierarchy of relations between correlation func- 
tions. We find that in the vacuum limit, where the system 
of equations describes few-particle scattering, the equa- 
tions for the dimer self energy and the splitting vertex 
are one-loop. Furthermore, the atom self energy is not 
renormalized. These three generalized couplings form a 
closed system of equations, decoupling from higher in- 
teraction vertices. We find the exact solution for these 
equations. These ingredients lead to the exact solution 
of the two-body problem, which manifests itself in the 
emergence of the non-perturbative Schrodinger equation 
for the bound state. 

In order to compute higher interaction vertices, we 
need to take higher loop diagrams into account, and the 
system of equations for these vertices is not closed. How- 
ever, we can establish a perturbative expansion of the 
equations in the limit J/\U\ 0. In [H], we establish 
the relation of the resulting effective theory to a spin- 1/2 
model in this limit. 



In this section, we analyze the quantum field theory 
derived above in the limits n = 0,2. This provides a 
useful starting point for the treatment of the many-body 
problem addressed in [14j . First we focus on the two- 
body problem, for which we present the exact solution 
within our framework. The Schrodinger equation for two- 
particle scattering is correctly reproduced. Clearly, for 
two particles, the physical three-body constraint cannot 



A. Dual Feshbach Model 

Before embarking the calculations, let us briefly dis- 
cuss the microscopic model emerging from our quantiza- 
tion prescription. We start with the microscopic action, 
obtained from the Hamiltonian Eq. (18) using Eq. (23 1. 



The complex fields t\ 2,^1,2 are now fluctuating classical 
variables and we may permute them at will, 



[tUr)t,,j{r)X,ir)X,{T) + 24_,(r)i2.,(r)4^^.(r)ii4T) 



(34) 



+\/2[4_,(T)<i,,(T)ii,,(r)+i2,.(r)4^,(r)4_^.(r)]X,(T) 



This action has the form of a Feshbach resonance model case) on the lattice, with non-local interaction parame- 
(see [29] for the fermion case, and [U [5] for the bosonic 
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ters: The dimer degree of freedom couples to the chem- 
ical potential with double strength, taking care of the 
double atom number in the dimer. The role of the "de- 
tuning" of the dimer state from the atoms is played by 
the onsite interaction. We note an interesting duality to 
the standard Feshbach model in the continuum for atom 
and dimer degrees of freedom obtained from a Hubbard- 
Stratonovich decoupling of an attractive (local) two-body 
interaction u: Such a procedure generates a detuning 
V ~ 1/u, while in our case the detuning v ^ U. Physi- 
cally, this means that the bound state formation in our 
lattice scenario is a weak coupling phenomenon, while 
being a strong coupling (resonant) effect in the contin- 
uum. Further note that the usual lattice construction 
using the single band approximation is delicate for sys- 
tems close to Feshbach resonances [3D] , and a realization 
of lattice Feshbach models in the resonant case is there- 
fore not a straightforward task, while it is realized in a 
natural way here. 

The physically most important coupling is the Fesh- 
bach or splitting vertex. It describes the formation of a 



dimer out of two atoms and its reverse. In contrast to 
the conventional continuum Feshbach model, the split- 
ting vertex in our model is bi-local, which in momentum 
space induces a form factor but does not lead to com- 
plications and still allows for an exact solution of the 
two-body problem. The term in the third line describes 
a non-local interaction between atoms and dimers. 

In the loop calculations we prefer to work in frequency 
and momentum space. With the definitions 



^'^'to.,„ tl,{T) = / e-'^--^,,, (35) 



x^ = (T,Xi), q = (w,q), 



2^ 



A=l 



the Fourier transformed action reads 



4,q(iw - M - 2eq) tl^q + t\^^{\UJ - 2fi + U) t2,q 

-V^ J S{qi -q2- 93)(eq2 + eq3)(4,qi*i,<72*i,?3 + h-c.) 

91,92,93 

-2 J S{qi - 92 + 93 - '74)[eqi-q, + eq2-q3]4,qi*l,924,93*2,94 
91, ...,94 

+ \/2 J (^(gi - 92 - 93 +94 - 95)[eqi-q2 +eqi-q3](4,gi^l,92^1,9.3^^94,9.'i + ^-C.) 
91, ■■■,95 

+ 2 J S{qi - 92 + 93 - 94)[eqi + eqj4.9i*l,92^^93,94 
91, ■■■,94 

- J ^iQl - 92 + 93 - 94 + 95 - 96)[eq2+q3-q4 + eq^+q5-q6]4,?i*l,92^^9..,94^^95 



(36) 



B. Dyson-Schwinger Equations 



Dyson-Schwinger equations (DSEs) [1^] are a direct 
consequence of the shift invariance of the functional in- 
tegral: 



= ^Jv{6x)^exp-S[x]+f5x 



(37) 



Switching to the effective action, i.e. requiring j = 
dV/Sx, the above equation turns into 



ST _ /SS\ 
6x \Sx/ 



]=5r/5x 



(38) 



This is the DSE for the one-point function. To reveal the 
structure of the DSEs for higher A^-point functions, we 
consider a general classical action with M vertices. We 
write the classical action in a vertex expansion. 



M ^ 

I ViSx) f - ^ + jVeicp-Six] + fSx. ^K] = S[x] + J2 wS^J!l.c.JXc.i ■ - • ^Xa«- (39) 

J ^ ox jv=i ■ 
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Here Ofj is a multi-index collecting field type as well as 
space and timelike (or momentum and frequency) indices. 
S and S^-^^ still depend on the classical field x- Plug- 
ging the vertex expansion into Eq. ( 38 ) relates the field 



derivative of the effective action to IPI Green functions 
up to order M. We can turn the DSE into a manifestly 
closed equation, i.e. an equation which is expressed solely 
in terms of the effective action and its functional deriva- 
tives. It reads 



ST 



c(l) I ^ q\-^) 



1 



(3) 



2! 



(40) 



M 



' (N - IV "i---ajv 



N-1 



s 



Gn 



For TV = 4 the derivative operator in the squared brack- 
ets is just the unit matrix. The full propagator is denoted 
by G and we have the relation Gap — (r^^^)^^- The full 
propagator as well as the classical vertices are functions 
of the classical field, S'(^) = S'(^)[x],G = G[x], such 
that the DSE for the A^-point correlation function can be 
obtained by taking iV — 1 functional derivatives on Eq. 
(40 1. We observe that the one-point function 5r/6xi3 de- 
pends on correlation functions up to order M — 1. Thus, 
the DSE for the TV-point function features vertices up 
to order M + TV — 2. Furthermore the self-consistency 
equations for the correlation functions for a theory with 
classical vertices up to order M features M — 2-loop di- 
agrams, since the order M vertex has M — I internal 
lines. 



C. Exact Solution of the Two-Body Problem 

The scattering problem is described by two coupled in- 
tegral equations for the exact dimer self-energy and for 
the exact Feshbach vertex, cf. Fig. [2] and App. |^ for 
the derivation. Inserting the frequency and momentum 
configurations appropriate for two-body scattering as de- 
picted in Fig. [2] and integrating out the frequencies, we 
obtain 

^(0)-l 



G^\E,k) = G'^^>-\E,k) + i:{E,k), 



E(i?,k) 

rk(p) 



1 f rk(q)r|f)(q) 



72 7 E + T^^\ci) ' 

dU rk(q)(rf (p) + r(°)(q)) 



rL^P) 



i? + rf (q) 



(41) 



with the definitions 

(q) = -2^(eq + eq_k), E = V2{iuj - 2/i), (42) 

where uj is the Euclidean external frequency. The 
difference between the full (Gd(-E,k)) and the bare 
{G'^P{E,'k)) Green's function is the dimer self-energy. 



s(£;,k) 

(a) 

k-p 

p 



k-q 

q 

k-p k p k-q 



vk 



FIG. 2: Scattering equations in vacuum. The solid lines 
represent atom propagators. Full vertices are signalled with 
heavy blobs, (a) Equation for the dimer self-energy. The full 
Feshbach vertex is needed for its solution, (b) Renormaliza- 
tion of the Feshbach vertex. External momentum configura- 
tions are chosen as needed in (a). 



!](£:, k) = G'^\£;,k)-G|j°^ ^{E,k). Note carefully the 
appearance of both external (p) and loop (q) momenta 
in the equation for the full Feshbach vertex. The second 
equation can be solved independently of the first one. 
This can be done by choosing the following ansatz for 
the full vertex, 

rk(q ) - T^^^ (q)7^°^ [E, k) + 7(1) [E, k), (43) 

where the unknown 7^"-' is dimensionless while 7*-^-' car- 
ries dimension of energy. The two unknown functions 
depend only on the external center-of-mass momentum 
k and the energy variable E\ the dependence on the rela- 
tive momentum q only appears in the coefficient of 7^°^ 
Comparing coefficients this ansatz yields the following 
system of coupled equations for the unknowns 



(0) 

(1) 



1 + 7(0' (1 -£;/) + 7(1)/, 



(44) 



1) 



7(")i?](i -£;/), 



where we use the abbreviation 

d'^q 1 



I{E,k)= J 



and the simplifications 



d^q T^°\q) 
d^q FL°)(q)Ff (q) 



(45) 



(46) 



i^^r E + T^°\q) 
The solution of the above equations is 

7(0) = 1, jW = E-r' 



E{EI - 1). 



(47) 



rk(q) = (£; + F(")(q))-/-i 
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The exact self-energy is given by 

!](£;, k) = -^(/-i-i?) = -(ic.-2^) 



(48) 



d'^q 1 1-1 

(2^)'i-2(eq + eq_k)+it^-2/iJ ' 



such that the equation for the full inverse dimer Green's 



function becomes, with G^*'^ ^(£',k) 



(2^^2(^ 



1 

fq-k) 



IW 



2/x 



Flu / 













(49) 



The presence of a bound state is signalled by a pole in the 
dimer Green's function at zero center-of-mass momentum 
and zero external frequency, G^^(a; = 0;/x,k = 0) =0. 
The chemical potential /i in the physical vacuum can 
be interpreted as the binding energy |28j after an ap- 
propriate decomposition which ensures that the atom at 
rest has no kinetic energy, fi = fit, — Jz. This defini- 
tion separates true kinetic from true potential (binding) 
energy, and finite momentum excitations have positive 
energy, Se^ = JJ2\i^~ cosqe^). Then, we have for 
the binding energy Eb = 2/x;,, and we note that the 
atoms are gapped out with half the binding energy, since 
their Green's function involves —fj-b- Thus, the molec- 
ular degrees of freedom are the lowest excitations since 
in contrast to the atoms they are massless (pole condi- 
tion). For vanishing binding energy we have fib = 0, such 
that the situation is reversed: The atoms are the gapless 
excitations, while the molecules are gapped. Introduc- 
ing dimensionless and dimensionally invariant variables 
U = U/{Jz), Eb — Eb/{Jz) the pole condition leads to 



1 

1^ 



d'^q 1 

-Eb + 2/dj:^{l - cos qe^) ' 



(50) 



This is precisely the Schrodinger equation for the dimer 
bound state. We discuss the formation of the bound 
state, since it will be interesting to confront these well- 
known results to the formation of di-hole bound states at 
n ~ 2 which shows different properties. The bound state 
forms at the critical Ud in d dimensions where Eb — 0. In 
three dimensions, the integral evaluates to a finite value, 
while in two dimensions a logarithmic infrared divergence 
pushes U2 to zero: 



U2 = 0. 



(51) 



This has to be compared to the mean field result at zero 
density, U^f — —2, obtained from the classical inverse 
dimer Green's function G'^°^^^(w = 0;^,k = 0). We 
observe a substantial downshift in the critical interac- 
tion strength. Close to the onset of the bound state the 
binding energy starts quadratically {d = 3) resp. expo- 
nentially {d = 2), due to the square root nonanalyticity 



FIG. 3: Dimensionless binding energy of dimers (n = 0) and 
di-holes (n = 2) as a function of the dimensionless interaction 
strength. The upper lines (black online) denote the mean 
field results, the lower/middle curves (red/blue online) are the 
exact binding energies in dimensions 2 and 3, respectively. 



resp. logarithmic divergence of the fluctuation integral: 

'\U\-\U3h^ 



d = 3 : Eb^- 
d = 2 : Eb^- 



CTU3 



(52) 



Sb,mf-C^mt-|C/|-2-|[/|, 

with numbers a ~ 0.42, A « 5.50 determined numeri- 
cally. We have added the linear mean field result. This 
dimensionally invariant behavior is approached for large 
negative couplings, cf. Fig. [3j 

The fact that the Schrodinger equation is repro- 
duced by our nonperturbative calculation is an impor- 
tant benchmark for our theory. Note that any deviating 
implementation of the constraint, such as an expansion 



of the square root Eq. ( 13 1 to leading order, would gener- 
ate incorrect prefactors in the scattering equations, and 
the Schrodinger equation could not be reproduced. 



D. Perturbative Limit: Effective Two-Body Hard 
Core Dimer Gas 

While the perturbative limit for the two-body prob- 
lem is straightforwardly obtained to any order from the 
exact solution above, for higher order vertices, nonper- 
turbative calculations are hard - the computation of the 
full dimer-dimer scattering vertex would require the com- 
plete solution of the four-body problem. However, in the 
limit J/U ^ 1 a systematic perturbative expansion of 
our set of Dyson-Schwinger Equations is available. For 
example, at order only diagrams with at most two ver- 
tices have to be taken into account. We find that a loop 
expansion to one-loop order is insufficient even at order 
J^. We may understand that qualitatively from the fact 
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that on short ranges, we are deaHng with a full quantum 
mechanical problem with strong fluctuations - the occu- 
pation of a site is either or 1. Thus, a loop expansion 
(in orders of h) cannot be expected to be reliable. 

In this section, we calculate the effective Hamiltonian 
for dimers. There dimers will obviously have a two-body 
hardcore constraint, as they are made up of two atoms 
each. The relation to a spin 1/2 model is made in [H] . 
There, we also show that a symmetry enhancement from 
the conventional U{1) ~ SO (2) for bosons to SO {3) is 
taking place in the strongly correlated limit, and discuss 
its physical implications. 

We will first calculate the effective Hamiltonian up to 
second order perturbation theory. We then perform a 



partial calculation of the fourth order. These results are 
crucially needed when discussing the many-body phases 
in the strongly coupled regime in [T3], as well as for the 
discussion of the nature of the phase transition. 



1. Second Order 

We are interested in extracting the effective action for 
dimers in the perturbative limit. If restricting to second 
order in J, the effective theory can only contain nearest 
neighbour terms. The most general form compatible with 
the constraint and symmetry principles is given by 



S[t2] = J dr(^z4,(r)(a,-Md)t2,»(T)-t^4_,(T)X,(T)X,(T)t2j(T) + ^^n^^^ 



r 



(53) 



The X terms are introduced in order to satisfy the con- 
straint principle, and there cannot be an onsite dimer 
interaction for the same reason. Z is a wave function 
renormalization factor. It accounts for the energy de- 
pendence of the perturbative expansion in x — J/\U\. 
fid = 2^ — U + Afid < is the effective dimer chemical 
potential, t, v are the constrained hopping and interac- 
tion constants to be determined. 

Kinetic Terms - The desired information on the local 
and the constrained hopping term can be extracted from 
the solution of the two-body problem, Eq. ( 49 ) . For this 
purpose we take 



-00, 



^ (w; /I, k) = iw - 2^ + [/ - 7 



4J^ 



y^(l + coske^). 



(54) 

In this form it is apparent that the full inverse propaga- 
tor contains the microscopic one exactly, but is now sup- 
plemented by a qualitatively new hopping term. To be 
consistent at second order, in the denominator we have 
to insert /i — —\U\/2. After Fourier transformation of 
the quadratic part of the action with inverse Green func- 
tion we obtain Afid = 2,Pz/\U\,t = 2,Pz/\U\. The 
wave function renormalization factor Z is extracted from 
Eq. (54 1 from expanding in iw; there is no second or- 



der contribution at second order perturbation theory, i.e. 
Z = 1. 

Dimer-Dimer Interaction - At short ranges, one ex- 
pects a dimer density-density repulsion due to the re- 
duced decay and recombination possibilities of one dimer 
if there is another one sitting close by. The Dyson- 
Schwinger Equation governing the dimer-dimer scatter- 
ing in the perturbative regime is displayed in Fig. [4] 
The one-particle irreducible graphs give rise to a nearest- 



neighbour (second order in J) density-density repulsion, 
due to the presence of the constraint, manifesting itself 
via the five-point splitting vertex. As anticipated above, 
it is interesting to note that the equation is two-loop even 
in the leading order perturbation theory. For the deriva- 
tion of the symmetry factors and the explicit calculation 
we refer to App. |A] Here we indicate the second order 
result of the momentum space calculation and discuss it: 



2 



2J^ 



f^ki -k2 + -k4 + £ki + fka + ^kg + ^ki ) -(55) 



The two qualitatively different momentum dependences 
correspond to different interaction processes in position 
space. After Fourier transform wc find the contribution 
to the effective action 



2J2 

W\ 



drY^ ["2,i('r)^2,i(T) + 4,j('^)*2,i(T)(n2,i -t- n2j) ■ 

(56) 



The terms may be interpreted as follows. The first term 
is a true dimer-dimer interaction describing the exchange 
of dimers on adjacent sites and is repulsive. The second 
one is the explicit manifestation of the constraint being 
inherited by the effective theory of dimers - they precisely 
contribute the terms linear in n2 which are contained in 
in Eq. 



(53) 



XiXj 

In summary, the effective couplings to second order 
read 



r. . . V 2.P 



(57) 



Both effective hopping and interaction are not present in 
the mean field approximation. Conceptually and prac- 
tically, they are however of high importance. The first 
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FIG. 4: The dimer-dimer interaction vertex to second order 
in the hopping J. 

one makes the dimers true physical, i.e. spatially prop- 
agating degrees of freedom, while the second one, with 
the positive sign, is very important for the many-body 
and long- wavelength calculations carried out in [14^, as 
it stabilizes the thermodynamic potential for the dimer 
superfluid. It also gives rise to the stiffness of the su- 
perfluid. Since the second order contributions are due to 
fluctuations on a single link of nearest neighbours, the 
results are dimension independent. 

2. Fourth Order 

At fourth order, we obtain not only nearest neighbour 
(nn) contributions to the constrained hopping and inter- 
action, but also next-to-nearest neighbour (nnn) contri- 
butions. However, due to the reduced number of path- 
ways connecting these more distant sites, the coefficients 
are substantially smaller then for nn terms. In the fol- 
lowing, we thus concentrate on the latter. In particular, 
for later purposes we will be interested in the deviation 
of the ratio of interaction vs. kinetic energy A = v/{2t) 
from the second order result A = 1. 

A brute force diagrammatic fourth order calculation is 
rather complex. Here we present a way to perform the 
fourth order calculation of the nn coeffcients for the in- 
teraction from a combination of geometric and diagram- 
matic arguments: First we argue that based on the ge- 
ometry of the contributing processes, the repulsive part 
of the fourth order contribution must equal the fourth 
order hopping contribution. The latter, in turn, can be 
calculated straightforwardly from the exact solution of 
the two-body problem. This argument can be applied 
for processes of arbitrary intersite distance, and we will 
present it in its general form. Then we refocus on near- 
est neighbours and identify processes contributing to the 
interaction which have no analog in the hopping process. 
We calculate the corresponding reduced set of diagrams 
explicitly. They yield an attractive contribution to the 
interaction strength such that A < 1. 

We begin with the geometric argument and consider 
the hopping first. If we are interested in the energetic 



contribution of hopping processes at a fixed order in per- 
turbation theory in J /\U\ and at a fixed distance, we can 
find this contribution in principle by drawing all possible 
pathways with a fixed number of hops which connect ini- 
tial and final site. We call these paths connecting initial 
and final state via hopping the hopping paths. Next we 
consider the interaction coefficient. The physical origin 
of this interaction emerges from the constraint: if there is 
a dimer at site j, then another particle cannot hop on this 
site. This can be calculated as the energy difference that 
emerges when comparing the number of paths through 
which a dimer at site i can decay and come back to this 
site without a dimer sitting at j, and the corresponding 
number when there is a dimer on j. The energy differ- 
ence obviously is always positive, because the constraint 
always excludes a stet of paths, leading to a repulsive in- 
teraction. Thus, the energy contribution can be obtained 
by just counting the number of paths where at least one 
of the travelling particles hits the site j. We call these 
paths interaction paths. 

Now we observe that any hopping path can be trans- 
formed into an interaction path by reversing the direc- 
tion of the arrows of one of the travelling particles on 
the shortest path which connects initial and final site 
(there may be several of these shortest paths) for any 
overall length of the paths, i.e. at any order of perturba- 
tion theory. Therefore, the number of interaction paths 
is larger or equal compared to the number of hopping 
paths. The reverse is also true. Thus the number of in- 
teraction and hopping paths is equal, and therefore the 
so-obtained contribution to the interaction must equal 
the fourth order hopping contribution. 

However, in general there are pathways contributing to 
the nn interaction which have no analog in the hopping 
pathways ^. These are processes in which none of the 
dimers is static, and we have to calculate them diagram- 
matically. For nn, the corresponding diagrams are pro- 
vided in Fig. [sj We discuss the corresponding processes 
(the explicit calculations are performed in App. |A|: In 
the left diagrams, the process starts with the decay of 
a dimer into atoms, where the sum over nearest neigh- 
bours indicates the z = 2d possibilities resulting from 
the unconstrained splitting vertex. The splitting is then 
followed by a double swap of atom and dimer. Finally, 
the atom on k recombines with the one on the target 
site i or j into a dimer. The position indices of these 
processes are fixed by the nn range of the couplings, and 
no further summation occurs. Performing the frequency 
integrations for zero external frequencies, the left side 
of Fig. [5] evaluates to —IQzJx^, with x = J/\U\, thus 
providing an attraction piece to the dimer-dimer interac- 
tion. The three-loop diagrams on the left take care of the 
constraint: one decay possibility described by the uncon- 



At second order, no other interaction processes can occur. This 
is the reason for A = 1 at second order. 
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FIG. 5: Attractive contribution to the nearest neighbour 
interaction at fourth order. The three-loop graphs implement 
the constraint, as discussed in the text. 



strained splitting vertex is not allowed, since the other 
dimer is located there. Indeed, the diagrams evaluate to 
+16 Jx^, such that the net result of the processes in Fig. 
El is 



^16(z- l)Jx^ 



(58) 



for real frequencies at a finite molecule momentum k, 
G^^(ia;(k);/ib = Eb/2,k) = 0, we find ^ 



4J 



(62) 



The wave function contributes a fourth order term to the 
dispersion relation. In order to see the effect on the ra- 
tio of interaction vs. kinetic energy, we absorb the wave 
function renormalization factor into a redefinition of the 
field t2 — > V~Zt2 and into the remaining couplings, e.g. 
t ^ i = t/Z,v ~^ V — v/Z"^, where the redefinition of 
the couplings is performed in such a way as to keep the 
full effective action invariant under the transformation. 
Since Z > 1 we conclude that A — v/{2t) < A, such that 
the ratio of interaction vs. kinetic energy is effectively 
decreased additionally. We finally note that the effects 
of the wave function renormalization correspond to an 
energy dependence of an effective Hamiltonian obtained 
in a Brillouin-Wigner perturbation theory for a micro- 
scopic Hamiltonian. It is well known that these effects 
occur at fourth order perturbation theory only. 



E. Relation of n = and n 
Mapping 



2: Particle-Hole 



This is an attractive contribution to the nn interaction 
in any dimension. 

We are now in the position to provide the full fourth 
order contribution. The repulsive part to the interaction 
is given by the fourth order contribution to the nn hop- 
ping, which we can obtain from the expansion of Eq. ( 49 ) 
up to fourth order with the result 



= t = 2J[x + {U{z - I) - 2)x^]. 



(59) 



We thus obtain the final result for the nn hopping and 
interaction terms. 



t = 2 J(a; + 2(6(z - 1) - l)^^), 

V _ Va+ Vr 

2 ~ 2 



(60) 



2J(a; + 2(2(z- 1) - l)^^). 



A = 



V _ 1 -h 2(2(z - 1) - l)a;2 
~2t ~ l + 2(6(z-l)-l)a;2 



1 



(z - \)x^ < 1. 



Finally, we discuss an additional effect at fourth order 
which is associated to the wave function renormalization. 
The leading term, relevant at fourth order perturbation 
theory, can be extracted from the frequency expansion of 
Eq. (l54|) at k = 0, yielding 



Z = 1 + 4zx^ 



(61) 



The correction to the dimensionless quantity Z is 0{x'^) 
and has effects at fourth order only. For example, 
if we are interested in the dispersion relation, ob- 
tained from the pole condition on the Green's function 



At density n — 2, the physics is expected to be similar 
to the case n = 0. Similar to the latter, the effect of spon- 
taneous symmetry breaking is absent due to the complete 
filling. Moreover, due to the constraint, there is only a 
single microscopic configuration, in full analogy to the 
completely empty lattice for n = 0. Importantly, there- 
fore the n = 2 state in the constrained model must not be 
viewed as a Mott insulator phase, such as e.g. n — I and 
strongly positive U, but rather as a constraint induced 
band insulator: For n = 1 the motion of the particles 
simply costs energy (providing for the gap in the Mott 
spectrum), while for n = 2 this is impossible in principle. 
Therefore, we may expect another "vacuum problem" . 
The low lying excitations on top of this vacuum will be 
holes and doublets of holes or di-holes instead of atoms 
and dimers, and the complete amplitude resides in the ^2 
mode. From this reasoning, we now replace 



t2.r ^ 1 



X'. 



(63) 



Alternatively, we could have adopted the formal point of 
view, applying the procedure described in Sec. |IV C| for 
the angles x = ^ = t''- Here we want to stress the physical 
similarities between n = 0, 2, but of course, the formal 
procedure generates exactly the same result. 



^ The propagation of the dimers in vacuum is gapless, and we have 
adjusted the zero of energy at k = via appropriate choice of /i. 
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With the above replacement, the Hamiltonian takes 
the form 

ffkin = -JY,[^th^^^'A.l+thtoAJl,J (64) 
Hpot — +1-1. 2tJ jto,i + tl.itld ~ ^0,i^O,i + 

+(-2/i + C/)Af^ 

where M is the number of lattice sites in each lattice 
direction. The potential energy term can be written as 

i i 

^i' = -^i + U. 

Working with the action in the following we may permute 
the operators at will. In this case, the mapping of the 
action for zero density to the one at n = 2 is given, up 
to constants (the mean field energies) by the following 
replacements 

t2 to, fi^ U ^U, 5spiit = V2J V2J, 

Jhop = J — > 2 J, ^exchange = 2J — > J. (66) 

While the coupling strength of the splitting vertex re- 
mains invariant under the transfo rmation, the roles of 
if^i") and iiT^^i) (^.f. Eqs. ( |l2|l8| ) are exchanged, and 
so is the role of the corresponding coupling constants. 
Based on Eq. ( 66 ) we observe that the scattering equa- 



tion for the holes can be obtained from simple replace- 



ments in the scattering equations (41 1. The equation for 
the hole splitting vertex, 



rk(p) 



rL°Hp) 



A rk(q)2(rf (p) + r;;^^(q)) 



(2^) 



i? + 2rf (q) 



(67) 



has the identical form to the one for dimers if we redefine 
the energy variable 

/2 

E ^ V2{iLu -2^1)^ E' [iu - 2/i') (68) 



and therefore has the solution (47 1 with the replace- 



(o)-i . 



ment E ^ E' . The bare inverse di-hole propagator is 
G|^°^"\£;',k) = iw + (-2^' + [/) =itj-(-2/i + C/), and 
the di-hole self energy E/j = G^^ — G 

I],(i?',k) = ll^[I-\E')-E') 



is found to be 
(69) 



iiuj - 2^') 



Bound state formation - The generalized "Schrodinger 
equation" , governing the formation of a bound state of 
holes, is given by the pole condition 



G^^(w 0;/i',k = 0) = 0. 



(70) 



With definitions analogous to the zero density limit (/i' — 
/i'j — 2Jz, E'j^ — 2fi'i^; this choice of /ij, ensures the zero 
of kinetic energy to appear at zero momentum), and in 
dimensionless variables, we find the explicit form 



1 



A\U\ - l2 + iEl 



d'^c 



{2TrY 



4/rfEA(l - cosqeA) 



(71) 



The bound state forms at the critical dimensionless inter- 
action strength defined by i?^ = 0. Due to the different 
structure of the scattering equations at n = 2, we expect 
different results from the n = case. Indeed, we find the 
dimensionless critical interaction strengths in three and 
two dimensions 
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-3.66, 



(72) 



An interesting feature occurs in d = 2: Despite the log- 
arithmic infrared divergence of the fluctuation integral, 
the di-hole bound state formation occurs at a finite at- 
tractive interaction strength. For the dependence of the 
binding energy on the interaction strength close to the 
threshold we find 



d = 3 : EL 



d = 2 : El 



m\-m) 

a(|^|-3) 



(73) 



-A^ exp( 



1(71-3 



), 



with (T, A given below Eq. (52 1. In contrast, the mean 
field binding energy reads Ei, = —\U\ +4. The exact 



results Eq. ( 72 ) for the onset of the di-hole bound state 
can be compared to the mean field answer Umi — —4. 
The shifts are substantially smaller than in the n = 
limit. Furthermore, from Fig. [3]we observe that the non- 
perturbative fluctuation dominated quadratic and expo- 
nential regimes are smaller than in the low density limit 
described by ( 52 1 . We conclude that fluctuations effects 
are weaker in the maximum density region n « 2 than in 
the low density counterpart n — 0. 

The finite value in Eq. ( 72 1 in two dimensions is sur- 
prising, and some comments are in order. An under- 
standing is obtained from the fact that the coupling con- 
stants are different in the band insulator n = 2 than 
in the vacuum n = 0, and so the coefficients in the 
effective propagators are different. This, in particu- 
lar, implies a shift in the effective interaction strength 
|C/| -> 4|f7|^2 + 3£'^; note that the integrals are identi- 



cal in Eq. ( 50 1 and Eq. ( 71 1 up to a coefficient 4 instead 
of 2 in front of the kinetic term in the integral - this 
refiects an enhanced mobility of the holes compared to 
the atoms. Therefore, we also have the log divergence in 
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d = 2 and n = 2 for i?^ — 0. However, this divergence 
now causes |[/| — 3 instead of |J7| — > for its compensa- 
tion. We emphasize that it is not a matter of definition 
of the interaction strength: [/ is a microscopic quantity, 
and the shift to it in Eq. ( 71 1 is due to the effects of 
the n = 2 band insulator. The finite threshold, thus, is 
an observable prediction of our theory. In addition, we 
would like to stress that there is no strong reason to ex- 
pect identical critical interaction strengths at n = and 
n = 2, as there is, in this regime of moderate interaction 
strengths, no particle-hole symmetry suggesting such be- 
havior. 

Effective di-hole theory - We have already noted that 
the mean field result for the binding energy is approached 
in the limit U — > — oo. The leading correction is given by 



1^1 V 



(1 + coske> 



(74) 



This precisely coincides with the leading perturbative 
self-energy correction for the dimers (cf. Eq. (54) 
at w = 0, — 2/Lt = \U\). Furthermore, the di-hole-di- 
hole interaction is identical to the zero density case at 
0{J'^/\U\), since the splitting and the five-point vertices 
entering the dimer-dimer scattering vertex (cf. Fig. |4]) 
have the same value, and only the atomic onsite propaga- 
tor enters at this order. Thus, to leading order we recover 
precisely the effective hardcore theory given by Eq. ( 53 ) 



with t2 (dimers) replaced by to (di-holes) but identical 
interaction constants, though in general the original mi- 
croscopic theories at n = 0, 2 feature different sets of 
interaction constants. The simple relation of the theory 
at n = and n — 2 must be expected: When the atom- 
atom attraction U is the largest scale in the problem, 
the integration of the high energy atom degrees of free- 
dom in perturbation theory can be performed prior to 
the inclusion of the effects of, e.g. a finite density in the 
system. Due to the decoupling from the atomic degrees 
of freedom, the theories in the perturbative regime but 
at different densities must then be directly related to a 
reference density, e.g. n = 0, by a rotation defined by 
Eq. ( 29 ) . The mapping of the theories for n = and 
n = 2 may be understood as performing the rotation 
with 9 = X — This is different from the case where 
the atoms do not decouple completely, such that the den- 
sity might have implicit fluctuation effects adding to the 
explicit effect of the rotation ( 29 ) on the dimer degrees 
of freedom. This hints at enhanced symmetry proper- 
ties in the perturbative limit, which are discussed in the 
companion paper [14 . 

Finally, we note that the low lying di-hole excita- 
tions on top of the n — 2 band insulator state disperses 
quadratically, as expected for nonrelativistic excitations. 
Once the density is lowered away from n — 2 and sponta- 
neous symmetry breaking sets in, there will be a linearly 
dispersing Goldstone mode |14j . 



VI. CONCLUSION 

In this paper, we have developed a method which al- 
lows to exactly map a bosonic lattice model with three- 
body onsite constraint to a theory for two unconstrained 
bosonic degrees of freedom with conventional polyno- 
mial interactions. The Gutzwiller mean field theory is 
recovered as zero order contribution to the thermody- 
namic potential. The quadratic fiuctuations around the 
mean field solution reproduce precisely spin wave the- 
ory. However, due to the exact nature of the mapping 
we can also systematically access the non-linearities. A 
convenient framework for our analysis is found to be the 
quantum effective action, the generating functional of the 
one-particle irreducible correlation functions. We estab- 
lish that the usual symmetry principles are supplemented 
with a new "constraint principle" , which depends on 
scale, being important at short distances but irrelevant at 
long wavelength. This setting allows to address fluctua- 
tions on various length scales in a unified framework. The 
application of the formalism to the full many-body prob- 
lem is performed in the related work [I4j , where we iden- 
tify various quantitative and qualitative aspects which 
are uniquely tied to the presence of interactions. Here, 
in order to demonstrate the validity of the formalism, 
and to prepare for the many-body analysis in |14) , we in- 
vestigate the scattering properties of few particles in the 
limit of vanishing density. Furthermore, we address the 
complementary problem in the limit of maximum filling, 
where the low lying excitations are holes and di-holes, 
and calculate the effective theory for hardcore dimers in 
the strong coupling limit. 

We believe that the formalism developed here has 
strong potential applications beyond the particular 
model analyzed here and in the related Ref. |T3]. We 
conclude by giving an outlook on problems that may be 
addressed within this framework. 

Exactly constrained theories - Spin models with spin 
S. The onsite Hilbert spaces can always be formulated in 
terms of 25* + 1 states. Here, we truncate to three states 
on each site {S = 1), but a generalization is straight- 
forward. Our construction maps such problems to 25" 
coupled bosonic degrees of freedom, and can thus be ef- 
ficient for small total spin S. More specifically, our for- 
malism can be particularly useful for XXZ models in 
strong external fields and a perturbative anisotropy. The 
mean field plus spin wave approximation is always incor- 
porated in our setting, but it also offers the opportunity 
to systematically assess the interaction effects. We also 
stress that inhomogeneous ground states, such as e.g. a 
charge density wave or antiferromagnet, are straightfor- 
wardly implemented in the formalism by applying site 
dependent rotations of the type (29 1. Furthermore, it 



will be interesting to investigate if the limit of infinite 
spatial dimension d ^ oo leads to a Gaussian fixed point 
of the theory, where the quadratic spin wave theory be- 
comes exact. This would add the possibility of a 1/d 
expansion to the conventional 1/5 (number of spin com- 
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ponents) and 1/A'' (number of field components in the 
corresponding nonlinear sigma model) expansion used to 
treat spin models. 

Models with an approximately realized constraint - A 
major challenge in quantum field theory is the treatment 
of strongly coupled systems where the non-linearity pro- 
vides a larger energy scale than the scales appearing in 
the quadratic part. Our theoretical framework is suited 
to address a certain class of such strongly coupled sys- 
tems, namely the one which is spanned by lattice models 
with strong local repulsion or attraction. It is not manda- 
tory that the onsite repulsion (two-, three- or few- body) 
is infinitely large. Often, only a few dominant low energy 
degrees of freedom are necessary to capture the essential 
physics, and these degrees of freedom are explicitly imple- 
mented in our approach. For example, in the Mott phase 
of the Bose-Hubbard model, number fluctuations around 
the Mott state with given n are small and a truncation 
to three Fock states n,n ± 1 is sensible ^Ij. Applying 
the formalism presented here, we are able to associate 
these particle and hole fluctuations with explicit bosonic 
degrees of freedom. In the resulting QFT, the nonper- 
turbative calculation of the phase boundary in the repul- 
sive Bose Hubbard model becomes feasible. Similarly, a 
strong but finite three-body repulsion 73 would lead to 
additional but massive bosonic degrees of freedom, which 
due to their mass ~ 73 can be taken into account pertur- 
batively. 

Fermions - The formalism is ideally suited to study 
attractive fermions on the lattice. In particular, recent 
experiments with 3-component fermions exhibit strong 
loss features [T^ and arc therefore important candidates 
for the observation of loss induced constrained models 
[7]. A modified scheme for the implementation of the 
constraint in one dimension has already been given in 
[7] , followed by the analysis of the model with bosoniza- 
tion techniques. Our scheme, which explicitly introduces 
a mean field describing qualitative features of the ground 
state, is promising in higher spatial dimensions, and its 
improvements compared to Gutzwiller mean field theory 
may be expected to be similar to Dynamic Mean Field 
Theory jMj . We note that an application to the fermionic 
repulsive Hubbard model is complicated by the choice of 
the qualitative features of the ground state being debat- 
able. In this treatment of the constraint along the 
lines of |35j may be preferrable. 
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Appendix A: Dyson-Schwinger Equations 

We calculate the DSE for the two- and three-point 
functions exactly, and the four-point function perturba- 
tively. In vacuum, a tremendeous simplification is pro- 
vided by the fact that density-type contributions vanish. 
In the loop language, this translates to checking whether 
there is a subdiagram in which the arrows form a closed 
tour. 

Formalism - The most important step is to find the 
correct symmetry factors for the loop contributions. For 
that purpose, it is sufficient to consider the field index 
a to parameterize the field type only; the frequency and 
momentum structure of the couplings can then be found 
from the corresponding conservation laws in a second 
step, reading off the couplings from the microscopic ac- 
tion ( 36 ) . Thus we consider 



Xa = (ii,il^2,4)- (Al) 
In this formalism, we find the following couplings: 



^ -5^^ = P^\^<^i&p2 + perm.] (A2) 
+P2['5q3^/J4 +perm.], 

where "perm." means to permute all greek indices while 
leaving the explicit numbers in their positions. This cou- 
pling must be inverted to give the propagators, with the 
result 



Gafi = —[5ai5i3i + perm.] + — [dazSfn + perm.]. (A3) 
Pi P2 



The higher couplings read 

^apj= [<5qi5/3i5',4 + Sa2Sp2S^3 + pcrm.] , (A4) 

^aPjSe= S^^''[^alSpiS^4.Ss4S^3 + (5Q2<5/32<573'554^e3 + perm.]. 

Pi,P2, s^'^K s^-^^ can be read off directly from the classical 



action ( 36 1 after proper symmetrization in momentum 



space. We can reduce the tensors by contracting with 
delta functions. For example, 

S'u4 = s''^^Sai6j3iS^4S''^^^ ^ s'-^^daiSpi[Sai6j3i +peTin.] 
^ 2s^^^SaiSai = 2s(3). (A5) 

where doubly occuring indices are summed over. The 
couplings are now fully symmetrized in field space, and 
indices may be permuted at will. 

Two- and three-point functions - We consider the 
DSEs for the two- and the three point functions for the 
inverse propagators and the splitting vertex. In the vac- 
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uum limit, these equations are one-loop: 





p(2) 


c(2) 


(2) 
12 




Plfi, 




r(3) 

^ 114 


_ c(3) 
— "^114 



Pi Pi 

-(3) 



1 1 



= 27(3)=2s(3)-4trs(4)AAy3)^ 

Pi Pi 

tr indicates the frequency and momentum integrations, 
while the discrete index contractions are carried out ex- 
plicitly. The corresponding integral equations for spe- 



cific external momentum configurations relevant for two- 
body scattering are given in Eq. (41 1 resp. Fig. [2] 

= 27^'^) holds under the assumption that the full 
Feshbach vertex has the same structure in field space as 
the classical one. The atom propagator is not renormal- 
ized in vacuum. Note that only the atom propagator 
enters the diagrams for the dimer self energy and the 
splitting vertex. 

Four-point functions - We further want to compute 
the fourth order interaction vertices. 
1) Second order, momentum space calculation - If we 
restrict to the order J^, we can limit ourselves to two- 
loop order. We find 



p(4) _ ,-,(4) rc.(5) ^ ^ -p(3) , c(5) (-< f< p(3) / a 7^ 

I c(^) p(3) _i_ p^^^ p(3) 1 

Jifr \<^'^^^ c n c r*-^' 4- i?^^' c c c r(^) 

31 L'-^^7MlA'2/i3 Wil'^l '^2/J2""P3l'3-^ aSviV^fi ^ *-'^(5miAI2P3 ^'I'^l '^2M2""P3l'3-^ (171/1 Iy2l'3 

I c(^) p(^) 1 

I 



Working at second order the full vertices are replaced 
by the classical ones. The dimer-dimer interaction then 
evaluates to 

rit'43 = 47(4) = _8tr 1 + 1 1^(5)] 

Pi Pi Pi Pi 

_8trs(5) 1^15(5). (A8) 

Pi P2,0 Pi 

The first two contributions are one-loop, the last one is 
two-loop. 



Now we insert the appropriate momentum structures. 
At order ,P the J dependence of the atom propagator 
is to be neglected and the dimer propagator equals the 
bare one, 

= iw - ^ = iw 4- |C/|/2, (A9) 
P2{'^) = P2,o('^) = iw — 2^ -I- [/ = icj. 

We then find, for zero external frequency, 



/" 1 1 

rit43 = +16J^( / —7-7: ^ I , X [iK-q + £q+q3-q4 + ^qi-q + gq+qi-q4][gq-q4 + ^q] + ^ q4:qi ^ q2}] 
I,,- pIu.)pX')p,{-1-u:') "^^'-^-- + '-^-^ + ^ q4}]^K-q + ^q.-q' + {Qi ^ Qal]) 



^qi + ^q2 ^ ^qs + ^q4 *^qi-q2 ^ ^qi-q4] ~ 



- 

The contribution to the effective action then reads 
r2 



2 J f 

AS" = y '5((7i - (72 +93 -g4))4('7i)*2('?2)4(93)i2('?4)[eqi +eq2 +eq3 +eq4 +eqi-q2 +eqi-q4] (AlO) 



91, ■■■,94 

r2 



2 J 

TTTf X! [4,»M*2j(T)[n2,j(r) -hn2j(r)] + n2,,(T)7i2j(T) 



2) Fourth order, real space nearest neighbours - Fi- nally, we compute the diagrams displayed in Fig. [5] As 
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discussed in the text, this is not the full contribution to 
the interaction vertex, which has a high diagrammatic 
complexity. (It could be obtained from Eq. (A7), to 
fourth order, but requires knowledge of F^^^ to second 
order.) However, knowing which diagrams we want to 
compute the procedure outlined above allows to identify 



the correct symmetry factors and contraction prescrip- 
tions. Since we are interested in a process on nn sites 
it is more efficient to do the calculation in frequency 
and real space, fixing the external legs at the nearest 
neighbour configuration, cf. Fig. [Sj We find 



-^(4)1.. . 

7 hr-A 



.2tr 11^(4) 1^(3) ll + 2tr \s^^^ 



Pi 



Pi P2 



Pi 



Pi 



= -16JS 
+I6J4 

= -f6(z- 1) 



duj 

did 
|[/|3 



duj' 1 



1 1 

Pi P2 
1 



1 1 

Pi P2 



(4)i__«(4). 
P2' 



= (3). 



Pi 



Pi 



(All) 



27r pi{u}) pI{~uj) Pi(w') p2{-uj - uj') 



duj' 
2^ 



duj'' 



2tt Pi{-uj - Lo') pi{uj") P2{(^ + uj' - uj") Pi{lu) P2{i^') Pi(w + a;') 



where the full vertices have been replaced by the bare 
counterparts andpi(u;) = iw — /i = \ijj+\U\/2,p2{uj) — ito, 
as appropriate at this order perturbation theory. The 
overall factor of two in the first line may be interpreted 
in terms of the two rows in Fig. [5| The factor of z m 



two-loop diagram results from the z possibilities for a 
dimer to decay via the bare splitting vertex. No further 
summations occur in real space. The three-loop diagram 
corrects for the fact that one of these possibilities is not 
allowed due to the presence of the other dimer. 
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